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ABSTRACT 

The photometric and morphological properties, as well as the star formation 
^0 I history, of the Draco dwarf spheroidal galaxy are analyzed on the basis of wide field, 

CCD photometry of the resolved stars covering about 1°^. 
' Draco is at a distance d = 80 it 7 kpc and has a metallicity, [Fe/H], of —1.8 it 0.2. No 

metallicity gradient is detected. The star surface density distribution can be fitted by 
a single exponential law of scale-length a = 5'.0 ±0'.l. The central surface magnitude 
is H'ly" = 24.4 lb 0.5 and the core radius Tc = 7'.5 ± 0^3 (equivalent to rc = 175 ± 7 
^ , pc). Within errors, the same scale-lengths are found for the density profiles along the 

semi-major axis and the semi- minor axis (re-scaled to semi-major axis units, using 
the ellipticity of the galaxy) of Draco. There are hence no evidences of a tidal-tail 
^ ' associated to Draco. The tidal radius of the galaxy is found to be rt — 42' (~ 1 kpc). 

d • The possibility that the large mass-to-light relation in Draco could be accounted 

for by a convenient spatial orientation is tested. An upper limit to Draco's size along 
the line of sight is ~ 14 kpc. This is too small to account for the velocity dispersion 
of Draco if it were due to projection effects only, and implies that other mechanisms 
(e.g., dark matter) are required. 

The stellar population of Draco is mainly old. Although some intermediate-age 
population is present in Draco, most of the star formation (up to 90%) took place 
before ~ 10 Gyr ago. No significant star formation activity is detected in the last ~ 2 
Gyr. Two methods {partial model and subgiant) have been used to investigate the star 
formation history of Drac, both producing results in good qualitative agreement. 

No difference is found between the scale lengths of the distributions of old (^ 9 
Gyr) and young (~ 2-3 Gyr) stars, indicating either that both populations were formed 
under the same kinematic conditions, or that any initial difference was afterwards 
erased. 
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Subject headings: galaxies: dwarf — galaxies: fundamental parameters — galaxies: 
individual (Draco) — galaxies: spheroidal — galaxies: stellar content — galaxies: 
structure 

1. Introduction 

In contrast with the early, simple idea that dwarf spheroidal (dSph) galaxies were systems 
essentially consisting of old, globular cluster-like, Population II stars (Baade 1963), they have 
been revealed as objects having undertaken complex star formation processes during their lives 
(Mighell & Rich 1996; Hurley-Keller, Mateo, & Nemec 1998; Stetson, Hcsser, 8z Smecker-Hane 
1998; Martfnez-Delgado, Gallart, Sz Aparicio 1999; Gallart et al. 1999b). Indications that dSphs 
could contain some intermediate-age or young stars were early found by Baade &: Swope (1961). 
Since then, other proofs have been accumulated, such as the presence of carbon stars (Aaronson 
& Mould 1980; Mould et al. 1982; Progel et al. 1982; Azzopardi, Lequeux, & Westerlund 1986) 
and bright AGB stars (Elston & Silva 1992; Frccdman 1992; Lee, Freedman, & Madore 1993; 
Davidge 1994; but see also Martmez-Delgado & Aparicio 1997). Stellar evolution theory became 
increasingly important in the interpretation of the color-magnitude diagrams (CMDs) when these 
reached the turn-offs of intermediate-age and old stars (see Mould &: Aaronson 1983 and Mighell 
1997 for the Carina dSph and Aparicio 1998 for a review). Comparison with theoretical isochrones 
allowed the determination of the range of ages present in the galaxies. DSph galaxies were studied 
with these techniques and the idea that all them, in one way or another, have a composite stellar 
population and that no two dSphs could be identified to have similar star formation histories 
(SFHs) became popular. In later years, the most powerful way to obtain the SFH, namely the 
analysis of the CMD with synthetic CMDs computed from a stellar evolution library (Aparicio 
2001), has confirmed the wide variety of scenarios and has revealed complex, frequently bursting 
SFHs in dSphs (Hurley-Keller et al. 1998; Mighell 1997; Martmez-Delgado et al. 1999; Gallart et 
al. 1999b; Hernandez, Gilmore, & Valls-Gabaud 2000; Aparicio 2001). 

Structurally, dSphs are usually well fitted by King models, although exponential fittings 
are also adequate (Irwin &; Hatzidimitriou 1995). They show relatively small scale lengths and 
have low central surface brightness (Gallagher & Wise 1994). They are the least-massive galaxies 
known, but still, their velocity dispersions imply mass-to-luminosity relations, T, as large as 
100 Tq. This is usually explained as dSphs containing large amounts of dark matter, but whether 
or not they are really bound systems (Gallagher &; Wise 1994) and whether or not the large 
values of T are due to peculiar orientation in space (Kroupa 1997; Klessen &; Kroupa 1998) is 
still subject of discussion. What is clear is that at least some dSphs, such as Sagittarius (Mateo, 
Olszewski, & Morrison 1998; Martmez-Delgado et al. 2001a) or Ursa Minor (Martfnez-Delgado 
et al. 2001b) show evidences of undergoing a process of tidal disruption by the Milky Way. This 
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would eventually be related to the galaxy formation predicted in cold dark matter scenarios, in 
which dwarf galaxies would be the first to form and would subsequently be merged into larger 
systems (White & Rees 1978). 

Several efforts have been devoted to the study of Draco. It was discovered by Wilson (1955) 
and was the subject of a careful photographic analysis leading to its first CMD by Baadc &: Swope 
(1961). They found several RR Lyrae variables and noticed that the horizontal branch (HB) was 
predominantly populated in the red part but also had a blue component. They were also the first 
to find anomalous Cepheids, later on identified in other dSph galaxies. Stetson (1979a) published 
new photographic photometry and estimated the distance modulus to be (m — M)o = 19.4 from 
the magnitude of the HB. Deeper CMDs were obtained by Stetson, VandenBergh, & McClure 
(1985), who reached the turn-off of the dominant, old population about 3.5 mag below the HB. 
Carney &: Seitzer (1986) found unambiguous evidence of a blue plume of stars above the main 
turn-off, as well as possible traces of an 8 Gyr-old stellar population. Recently, Grillmair et al. 
(1998) have presented WFPC2 HST photometry that allows the MS luminosity function down to 
3 mag below the oldest turn-off to be determined. 

The stellar density distribution of Draco was early studied by Hodge (1964) and, more 
recently, by Irwin &; Hatzidimitriou (1995) who report a tidal radius of 28' and by Lake (1990), 
Pryor Sz Kormendi (1990) and Piatek et al. (2001) who find evidences of a larger extension of the 
galaxy. 

Draco's metallicity has been found to be very low. Prom the color of the RGB, Carney 
& Seitzer (1986) estimated it to be [Fe/H] ~ —2 with a dispersion of 0.8 dex. This is in good 
agreement with spectroscopic results for several giant stars, which lead to —3 < [Fe/H] < —1.5 
(Stetson 1984; Shetrone, Bohe, & Stetson 1998); [Fe/H] = -1.9 ± 0.4 (Lehnert et al. 1992) or 
[Fe/H] = -2.00 ±0.21 (Shetrone, Cote, & Sargent 2001). This studies also clearly establish a high 
metallicity dispersion of the order of 1 dex. 

Since the early work by Baade &; Swope (1961), the SFH of Draco has been the subject of 
much discussion. The position and morphology of its RGB would indicate a very old age for 
the main population, similar to that of the globular clusters (Grillmair et al. 1998). But the 
predominantly red HB seems to be incompatible with this, at least if age is the second parameter 
(although recent results indicate that this assumption might be inadequate: Rosenberg et al. 1999; 
Behr et al. 2000). In such a case, Draco should be some 2 Gyr younger than both the globular 
clusters and the Ursa Minor dSph galaxy (Stetson et al. 1985; Carney & Seitzer 1986). On another 
hand, Carney &; Seitzer (1986) found some evidence of a younger, ~ 8 Gyr-old population and, if 
the blue plume is produced by younger MS stars rather than by blue stragglers, in favor of the 
hypothesis that the star formation in Draco should be extended up to a few Gyr ago, even if at 
a low rate (Carney & Seitzer 1986; Grillmair et al. 1998). This would also be compatible with 
the presence of the aforementioned anomalous Cepheids, which would belong to a few Gyr-old 
population (Norris & Zinn 1975; Demarque & Hirshfeld 1975; see also Gallart et al. 1999a and 
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references therein). 

In this paper we present an analysis of the star formation history, metahicity, and 
morphological properties of Draco, based on wide-field, ground-based photometry. We confirm 
previous results favoring a predominantly old stellar population and show evidence of a relatively 
short star formation burst 2-3 Gyr ago. The paper is organized as follows. In §2, the observations 
and the data reduction are presented. In §3, the CMD of Draco is presented and the metallicity 
and distance obtained. In §4, the morphological properties of Draco are discussed, including its 
overall structure, surface luminosity distribution and three-dimensional orientation. In §5, the star 
formation history is computed using two different methods. Finally, the main results of the paper 
are summarized in §6. 



2. Observations and data reduction 

Images of Draco were obtained in B and R Johnson-Cousins filters during several observing 
runs with the INT (2.5 m) at Roque de los Muchachos Observatory on the island of La Palma 
(Canary Islands, Spain). Table || gives the journal of observations. The WFC was used with four 
2048 X 4096 EEV chips. The scale is 0.33 "/pix which provides a total field of about 35 x 35 
arcmin^. Several fields were sampled, as shown in Fig. |^. Figure |2| displays an isopleth map of 
Draco, by Irwin & Hatzidimitriou (1995), showing also the area covered by our observations. Total 
integration times were 1920 s in B and 1600 s in i? for field A (the central field); 3600 s in i? and 
4200 s in R for field B, and 2820 s both in B and R for field C. 

Bias and flat-field corrections were made with IRAF. Then, DAOPHOT and ALLFRAME 
(Stetson 1994) were used to obtain the instrumental photometry of the resolved stars. 
Transformation into the standard photometric system requires several observations of standard 
star fields in such a way that standards are measured in all four chips. In practice, we observed 
three fields in the four chips and used them to calculate relative photometric transformations from 
each chip to the central one (chip 4). Then, a larger number of standard star fields were measured 
in chip 4 during the observing run to calculate both nightly atmospheric extinctions and the 
general transformation into the standard Johnson-Cousins system. In total, 280 measurements of 
30 standards were made during the observing run of June 1999. Summarizing, the transformations 
from each chip to chip 4 are: 



{bi - bi) = 0.36 - OM{B - R); cr = 0.01, (1) 

(r4-ri) = 0.45 + 0.04(5 -i?); (7 = 0.02, (2) 

(64-62) = 0.409 -0.003(5 -i?); cr = 0.007, (3) 

(r4-r2) = 0.462 + 0.102(5 -i?); = 0.003, (4) 

(64-63) = 0.438 -0.024(5 -ii); cr = 0.008, (5) 
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and 

(r4-r3) = 0.48 + 0.06(5 -i?); a = 0.03, (6) 

where subindices refer to chips, lower-case letters stand for instrumental magnitudes, and capital 
letters for Johnson-Cousins magnitudes. The a values are the dispersions of the fits at the centers 
of mass of the point distributions; hence they are the minimum internal zero-point errors. There is 
a large zero point (bigger than 0.4 mag in most cases) displacement between chip 4 and the other 
chips, which implies the lower sensitivity of the former. 

Transformations from chip 4 instrumental magnitudes measured at the top of the atmosphere 
into Johnson-Cousins magnitudes are given by: 

{B - 64) = 24.603 + 0.036(5 - i?); a = 0.006 (7) 

and 

(ii:-r-4) = 24.504 + 0.040(5 -i?); a = 0.010, (8) 

where, as before, lower-case letters refer to instrumental magnitudes and capital letters to 
Johnson-Cousins magnitudes. Finally, dispersions of the extinctions for each night are about 
a = 0.01 for both filters. 

This provides the photometric transformation for the 1999 observing run. For the remaining 
runs, transformations were done using as secondary standards the stars of the Draco fields 
overlapping with those of the 1999 run. In this way, the consistency of the photometric scales of 
the different fields is secured. 

Aperture corrections are necessary to transform the ALLFRAME magnitudes into large-radius 
aperture magnitudes. The latter is the system in which the magnitudes of standard stars have 
been measured. The transformation has been performed through aperture photometry on a set 
of ~ 200 isolated bright stars of field A after subtracting the remaining stars from the frames. 
The difference between these and the ALLFRAME magnitudes for the same stars provides the 
aperture correction for the corresponding frame. Dispersions in the aperture corrections are of the 
order of 0.01. Putting all the errors together, the total zero-point error of our photometry can be 
estimated to be about a = 0.02 — 0.03 for both bands, although it is up to cr = 0.04 for chip 3 in 
filter R. 

For each star, ALLFRAME provides with the dispersion in the PSF fitting. In general, these 
dispersions do not reproduce the external errors of the photometry, which stem mainly from 
stellar crowding and blending (see Aparicio, &: Gallart 1995 and Gallart, Aparicio, & Vflchez 
1996), but they provide an indication of the internal accuracy of the photometry. In this sense, 
our average ALLSTAR errors are 0.02, 0.06, and 0.15 for R = 21.1, 22.8, and 23.9, respectively, 
and for B = 22.5, 23.9, and 25.0, respectively. 
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Artificial-star tests were performed on field A in the usual way (Stetson & Harris 1988; 
see also Aparicio &: Gallart 1995) to check the observational effects and estimate completeness 
factors as a function of magnitude. In practice, 20 000 artificial stars were added to the B and 
R images of Draco. Coordinates of these stars were distributed in a grid covering the whole 
image. The step used to build this net was 7" both in x and y, avoiding overcrowding effects. 
Colors and magnitudes of artificial stars were obtained from the synthetic CMD (computed as 
described in Aparicio & Gallart 1995) of an old stellar population. In this way, they approximately 
mimic the real CMD of Draco, which is important for a realistic observational effect test. The 
resulting completeness factors are shown in Figure ^. They start to drop at "V" ~ 23.0 (with 
"V" = {B + R)/2; see §3.1). But about 5% of the stars are lost at any magnitude brighter than 
"V" ~ 23.0. These are mostly stars lying in the regions around very bright, saturated, field stars. 

3. The stellar content, metallicity and distance of Draco 

3.1. The color magnitude diagram 

Figure ^ shows the CMD of Draco. Stars within the inner 30' (the semi-major axis) have been 
plotted (see below). A pseudo-F magnitude, obtained as "V" = {B + -R)/2, is used in the vertical 
axis. This has the advantage of producing an almost horizontal HB. It shows a well populated, 
narrow red giant branch (RGB) with a relatively small slope as well as a blueward extended HB. 
These features are the trace of an old, low-metallicity stellar population with a small metallicity 
dispersion. The main sequence (MS) is visible at "V" > 23 and, very interestingly, an extended 
sequence of blue stars with "V" < 23, {B — R) < 0.6 is also present. This could be a blue-straggler 
(BS) sequence or the signature of a younger stellar population. 

Figure ^ shows CMDs for elliptical annuli of increasing semi-major axis. Note firstly that 
traces of Draco RGB and MS are still evident at 20' < a < 42'. Also noticeable is the small group 
of stars above the red part of the horizontal branch, visible in the CMD for a < 7/ farcm5 at 
"V" ~ 19.7, {B — R) = 1.0. This can be a red clump formed by stars burning He in their nuclei 
and would be the signature of a young-to-intermediate population. In fact, Baade & Swope (1961) 
and Zinn & Searle (1976) found anomalous Cepheids in this region. An explanation for the origin 
of these variables is indeed that they are young-to-intermediate age stars. If so, they should have 
an MS counterpart, perhaps like the blue extension discussed above. We will come again to this 
point in Section 4. 

These CMDs will be used later on as a reference for discussion of galactocentric gradients of 
stellar populations and metallicity. 
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3.2. The metallicity 

The metallicity of Draco has been investigated by Carney &: Seitzer (1986), who obtained 
[Fe/H] = —2 from the color of the RGB. From the spectroscopy of three giant stars, Shetrone et 
al. (1998) found that the metallicity should be in the range —3.0 < [Fe/H] < —1.5. 

The color of the RGB is frequently used to obtain an alternative estimate of the metallicity. 
Recently, Gallart et al. (2001) have shown that these values can be wrong. However we have 
applied the method to our data and derived an estimate of the metallicity in the understanding 
that it will be useful for comparison with results for other galaxies existing in the literature. 

To do so, we have used the relations provided by Saviane et al. (2000), based on the position of 
the RGB of several globular clusters. However, these relations are calculated in the [My, {V — I)] 
plane while our photometry is BR. Hence the Saviane et a/.'s relations must be transformed into 
the [M/j, [B — R)] plane to be used with our data. 

For this we have used the Padua (Bertelli et al. 1994) and Yale (Yi, Demarque, & Oemler 
1997) isochrones to obtain, for a set of metallicities, the values of the absolute Mr magnitude and 
the color indices corresponding to the RGB at Mj = —3. We denote these values by Mrj^_^, 
{V — I)/=_3 and {B — i?)/=_3, respectively. In this way, we obtain the curve in the [M/j, {B — R)] 
plane shown by a thick line in Figure ^. Besides Mjij^_^ and {B — i?)/=_3, the parameter changing 
along this curve is (V — I)i=-^. The intersection of the RGB fiducial line of Draco with the curve 
gives its iy — /)/=-3, whose value is 1.18. When this is introduced into the Saviane et al. (2000) 
relations, the metallicity is obtained to be [Fe/H] = —2.0 it 0.1 if the Zinn & West (1984; ZW from 
here on) metallicity scale is used, or [Fe/H] = —1.6 ±0.1 if the Carreta &l Gratton (1997; CG from 
here on) metallicity scale is used. 

Besides this, we have directly obtained an alternative calibration for the metallicity in the 
[Mr, {B — R)] plane. For this purpose, we have used the RGBs of three globular clusters: NGC 
1904 (Alcaino et al. 1994), NGC 4590 (Alcaino et al. 1990), and NGC 6723 (Alcaino et al. 1999). 
A summary of the metallicity (according to scales for the globular cluster metallicities by ZW 
and CG), the distance and the reddening of these clusters is given in Table ^. In accordance with 
Saviane et al. (2000), the cluster RGBs have been fitted by a hyperbolic equation of the form 
y = a + bx + c/{x — d). This fit has then been used to calculate the color index {B — R) for 
M/j=_2.5, which we will denote {B — R)ji=-2.5- Finally, the relation between the metallicity and 
{B — R)r=-2.5 has been determined for both the ZW and CG scales. For the second case, a linear 
fit of the form [Fe/H] = a{B — R)r=-2.5 + b is enough, but for the ZW scale, a quadratic fit of the 
form [Fe/H] = a{B — R)\^_2 5 + b{B — R)r=-2.5 + c is preferable. The coefficients of the fits are 
given in Table ^ 

The metallicities obtained are [Fe/H] = —1.9 it 0.1 if the ZW metallicity scale is used or 
[Fe/H] = -1.6 ± 0.1 if the CG scale is used. 

Considering all the values obtained here, we adopt [Fe/H] = —1.8 it 0.2 for the metallicity of 
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Draco. Figure |^ shows the RGBs of the three clusters used to calculate the [Fe/H] —{B — R)r=-2.5 
relation superimposed on the CMD of Draco, and demonstrates that, in fact, the metallicity of the 
galaxy should be between the NGC 4590 and NGC 1904 values. 

(B — R)r=-2.5 and the metallicity obtained therefrom can be used to test whether 
galactocentric metallicity gradients exist in Draco. CMDs for elliptical annuli of increasing 
semi-major axis, similar to those shown in Fig. |5|, have been used to this purpose. Third-degree 
polynomial fits to points in the interval —1.5 < R < —3.2 were used to obtain a fiducial RGB for 
each CMD. The (B — R)r=-2.5 values for each one are represented in Figure |[ The error bars 
show the internal dispersions of points about the polynomial fits. Inspection of this figure reveals 
that no metallicity gradient can be deduced to exist in Draco. 

A metallicity dispersion of about 1 dex in [Fe/H] has been clearly established in studies 
referred to the inner ~ 5' of Draco (Stetson 1984; Shetrone, Bolte, & Stetson 1998; Lehnert et al. 
1992; Shetrone, Cote, & Sargent 2001). Note that this is not in contradiction with the absence of 
an overall, large-scale metallicity gradient found here. 

3.3. The distance 

Previous estimates of Draco's distance have been made by Stetson (1979a), who obtained 
(i = 72 it 3 kpc using the position of the red HB, and by Nemec (1984), who obtained d = 84 it 12 
kpc using RR Lyr variability periods. Using our data, the distance to Draco can be estimated 
through the magnitude of the HB at the RR Lyrae instability strip. From an analytical fit to the 
blue part of the HB we have estimated that magnitude to be "V" = 20.2 it 0.1. A shift has to be 
applied to this value to obtain the actual V magnitude. Considering the BVR magnitudes in the 
RR Lyrae instability strip given by the Padua stellar evolution library, that shift is "V" — V = 0.06, 
hence V = 20.14. Finally, using the relation by Lee et al. (1999) for the luminosity- metallicity 
of the HB, My = 0.17[Fe/H]+0.82, the distance becomes (m - M)v = 19.6 ± 0.2, where the 
quoted error includes uncertainties in the "V" — V transformation, as well as the metallicity error. 
Reddening for Draco was investigated by Stetson (1979b), who found E(B — V)= 0.03 ± 0.01. It 
can also be calculated from the IR maps of dust by Schlegel, Finkbeiner, & Davis (1998), which 
give E{B — V) = 0.028, corresponding to an extinction of Ay = 0.09. Applying this correction to 
the distance modulus gives (m — M)o = 19.5 ± 0.2, or d = 80 ± 7 kpc. We will adopt this value. 

4. The morphological properties of Draco 

4.1. Density profile 

The radial stellar density distribution of Draco was first obtained from wide-field photographic 
photometry by Hodge (1964), who reported a tidal radius rt = 26', and more recently by Irwin & 
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Hatzidimitriou (1995), who derived rt = 28'. Lake (1990) and Pryor & Kormendi (1990) modeled 
the galaxy using Hodge's data and found rt = 55' and rt = 36', respectively. Very recently, Piatek 
et al. (2001) have carried out an extensive search of Draco stars beyond the Irwin & Hatzidimitriou 
(1995) tidal radius, finding evidence of the existence of such stars. 

Here we provide a new profile, now based on CCD photometry. Three stellar distributions 
have been obtained from star counts in elliptical annuli of increasing semi-major axis. The 
eccentricity of the elliptical annuli has been assumed to be e = 1 — h/a = 0.29 and the position 
angle 82°, after Irwin &; Hatzidimitriou (1995). The center of the galaxy given by these authors has 
also been adopted. The three density profiles have been computed as follows. For the first one, all 
the resolved stars have been used. For the other two, only stars in two regions of the CMD where 
the ratio of Draco to foreground stars is largest, have been considered: the MS and the RGB. To 
sample the MS, stars within the box defined by 0.6 < {B - R) < 1.0; 24.25 <" V" < 22.75 have 
been used. For the RGB, box 14, as defined in Figure has been used. 

The final star density distributions are shown in Figure ^. The corresponding "V" surface 
magnitude is given in the right-hand vertical scale, according to the magnitude scale obtained 
below. In all the cases, open circles show the raw distributions. The sky contamination by 
foreground and background objects must be subtracted from this to obtain the distribution of 
Draco stars. Background, extended objects are automatically removed by DAOPHOT, while 
background, stellar-shaped objects are indistinguishable from foreground stars. For this reason, 
from here on, we will refer to foreground contamination only, understanding that it include also 
the background stellar-shaped objects not rejected by DAOPHOT. 

The raw radial surface density distributions (Fig. open circles) change slope at a ~ 20' — 25', 
but continue decreasing until at least a ~ 42'. Let assume that the sky star density has been 
reached at this distance. If the average densities for a > 42' are subtracted from the raw radial 
distributions, the profiles shown by full dots in Figure ^ are obtained. A single law can account for 
the density profile for a > 5' in all the cases. On the one hand, the absence of a plateau for large 
galactocentric distances indicates that a tidal tail is not present in Draco. On the other hand, this 
indicates that the tidal radius (semi-major axis) of Draco is about rt = 42'. 

A way to further explore the presence of extra-tidal stars or of a tidal tail structure is to 
compare the density profiles in the direction of the major and minor axes. We have done so for 
the MS and the RGB stars, for which the Draco-to-foreground star ratio is greater. The results 
are shown in Figure The distributions obtained from star counting in full ellipses (the same 
shown in Fig. full dots) are compared with those for stars within two strips along the western 
semi-major axis and the southern semi-minor axis. Table ^ lists the scale- lengths in each case 
(column 2). For the semi-minor axis, these quantities have been scaled into units along the 
semi-major axis; i.e., they have been divided by b/a, so that they are directly comparable with 
the scale-lengths corresponding to the semi-major axis and full ellipses. No significant differences 
are found in either case, which leads to the conclusion that Draco is a very extended, low density 
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galaxy, but that no evidence of a tidal tail populated by stars stripped away from the galaxy are 
found. In any case, more data about the stellar population in fields around Draco's main body are 
required to better understand the galaxy's structural properties. 

Finally, the core radius (semi-major axis) is Vc = 7f5ib0!3. At the distance adopted for Draco, 
this corresponds to Vc = 175 ± 7 pc. 

4.2. Surface photometry 

It is difficult to obtain the surface brightness distribution of galaxies which, like Draco, have 
a very low stellar density and which are resolved down to stars of quite faint intrinsic luminosity. 
An indirect technique is to obtain the surface density distribution of resolved stars, which is 
afterwards transformed into a surface-brightness scale. In our case, this transformation has been 
done as follows. Firstly, a synthetic model for a stellar population as found for Draco in §5 has 
been computed. Secondly, the number of RGB stars has been counted in this model. To do this, 
stars in box 14 (see Figure ^) have been considered. Finally, comparison of this number with the 
total luminosity in B and R for the synthetic model provides a scale transformation from RGB 
star counts to magnitudes. The accuracy of the final scale of magnitudes cannot be expected 
to be high. In particular, it depends on the assumed IMF. But it is probably better than that 
obtained through direct measurement of brightness on a big field dominated by background noise 
and contamination by bright foreground stars. 

Applying this method, we estimate the central surface brightness in B, R, and 'V to be 
fj,B,Q = 25.0±0.5, fiRfi = 23.9±0.5 and H"v",o = 24.4ib0.5, where the quoted errors should account 
for the high degree of uncertainty. Irwin & Hatzidimitriou (1995) obtained nv,o = 25.5 ± 0.5, 
which is one magnitude fainter than ours. It is difficult to say which value is preferable. Irwin 
&: Hatzidimitriou have obtained theirs from a King model fit to photographic data. Our CCD 
photometry is deeper and is expected to be more accurate, but, as previously quoted, our estimate 
could be affected by a wrong choice of the IMF and, in general, by the fact that we are using a 
stellar population model. In our opinion, the relatively large disagreement between both values 
reveals nothing but the difficulty of obtaining accurate surface photometry of this kind of systems. 

The integrated magnitudes of Draco can be obtained by integration of the density profile and 
the central surface brightness abovely obtained. The results, together with all the morphological 
parameters obtained here and several other integrated quantities are summarized in Table ||. In 
Table ^, morphological parameters are compared with other authors' results. 
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4.3. Spatial distribution of stellar populations 

The spatial distribution of stellar populations of different ages gives information on the 
time-scale in which stars mix out in the galaxy and on kinematic and/or morphological evolution. 
Computing the detailed models required to perform an accurate analysis is beyond the scope 
of this paper. But valuable information can be obtained from the comparison of the radial 
distribution scale lengths of populations of different ages. The latter can be selected from the 



boxes defined in Fig. 13, since different boxes sample different ages (see details in §5). 



From synthetic models computed for the SFH study (see below, §5) it can be seen that boxes 

12 and 14 (Fig. [l^) contain stars with ages in the interval 9-15 Gyr, while boxes 2 to 6, 10 and 

13 contain younger stars and preferently, stars in the interval 2-3 Gyr (they sample the small 
recent starburst; see §5). We have calculated the radial distribution scale length for both groups 
of stars in the region 5' < a < 20'. We have obtained a = 5'.8 it 0'.2 for the older population and 
a = 6' lb 1' for the younger; i.e., there is no significant difference either between both or with 
respect to the whole resolved stellar population. Within the errors, this indicates that stars of 
different ages where formed with similar spatial distribution or that their movement in the galaxy 
has completely mixed them. 



4.4. On the three-dimensional orientation of Draco 

Dwarf spheroidal galaxies usually present high mass-to-luminosity ratios (T). Draco, in 
particular, has T = 84 Tq (Mateo 1998). Dark matter is usually invoked as an explanation 
for these large values in a virialized scenario. An alternative explanation has been proposed by 
Kroupa (1997) and Klessen & Kroupa (1998), who suggest that tidal disruption of these galaxies 
by the Milky Way potential would well produce very elongated systems whose main axis, together 
with the main axis of the velocity ellipsoid, could eventually be oriented close to the line of sight. 
In such a case, the high observed velocity dispersion would not need of large quantities of dark 
matter to be accounted for. 

To check this possibility we have adopted the hypothesis that Draco is configured according 
to model RSl-4 of Kroupa (1997), which is the one that most closely reproduces the T of Draco. 



The spatial distribution of stars in this model is shown in Figure [11]. We have then simulated 
the predicted distance dispersion of Draco stars to us in a synthetic CMD of a stellar population 
similar to that found for Draco (see §4). The CMDs before and after simulating the distance 
dispersion are shown in Figure |l^. Comparison of the right-hand CMD with the Draco CMD (Fig. 
^) clearly indicates the incompatibility of the RSl-4 model with the observations. 

On another hand, the extension of Draco along the line of sight can be estimated from the 
width of the HB. This width is potentially affected by several factors, such as observational effects, 
stellar evolution or binarity of composite stellar populations. From our artificial-star test (§2) we 
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have checked that widening of the HB by observational effects can be neglected, but considering 
every other possible effect in detail is complicated. However, an upper limit estimate can still be 
set by neglecting all the effects and assuming that the entire width observed in the red HB (Fig. 
^) is produced by distance dispersion. In such a case, the size of Draco along the line of sight 
would be ~ 14 kpc. This is quite a large value compared with Draco's size in the plane of sky but 
it should be stressed that it is an upper limit, and that it is still less than 1/3 the size given by 
the RSl-4 model. 

5. The Star Formation History 

The SFH of a galaxy can be derived in detail from a deep CMD through comparison with 
synthetic CMDs (Aparicio 2001). Details on the computation of synthetic CMDs and on the 
method used here to derive the SFH can be found in Gallart et al. (1999b), Aparicio &: Gallart 
(1995), and references therein. The accuracy and temporal resolution of the SFH depend on the 
quality and depth of the observed CMD and on the uncertainties in the stellar evolutionary library 
upon which the synthetic CMDs are computed. To obtain a reasonable resolution for old ages 
the CMD must sample stellar evolutionary phases sensitive to these ages, such as the HB or the 
lower MS. Moreover, the RC, the RGB, and the HB itself provide information on the chemical 
enrichment law (CEL). 

The CEL and the initial mass function (IMF), together with the star formation rate (SFR), 
are the functions contained in the more general SFH. A function containing information on the 
fraction and mass distribution of binary stars would also be relevant. For simplicity, we have 
neglected here the influence of binary stars and have assumed a fixed IMF: that given by Kroupa, 
Tout, & Gilmore (1993). The CEL can be constrained by the morphology and color of the RGB. 
As discussed above (see §3.2) Draco presents a low average metallicity and a low metallicity 
dispersion, [Fe/H] = -1.8 ± 0.2. This corresponds to 0.0002 <Z< 0.0003. However, the Padua 
stellar evolutionary models produce a best fit to Draco RGB for Z ~ 0.0005. This difference is due 
to the internal tuning of the stellar evolution models and is not important for the analysis of the 
metallicity. However, it is preferable using it to compute the synthetic CMDs necessary for the 
SFH study. In practice, we have assumed a CEL independent of time producing metallicity values 
for the stars randomly distributed in the interval 0.0004 < Z < 0.0006 for all ages. Provided that 
the stellar population does not have a large amount of young stars (younger than some 1 Gyr), as 
can be deduced from the extension of the MS of Draco, this CEL produces synthetic CMDs with 
an RGB always compatible with that of Draco. 

The only remaining function to complete the SFH of Draco is the SFR. Two different methods 
have been used to obtain it. The first method is that of the lAC-Padua group (see Aparicio 1998 
for a description and Gallart et al. 1999b for one of the latest applications), which we will refer to 
as the partial diagram or partial model method. The second method is that presented by Mighell 
(1997), which we will refer to as the subgiant method. We will discuss in turn the results obtained 
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with each of these. 

5.1. The partial diagram method 

This method was introduced by Aparicio, Gallart, & Bertelh (1997). To apply it, a few 
synthetic CMDs (partial models) are computed, each for a stellar population with ages in a 
narrow interval (typically, a few partial models are used each embracing an age interval ranging 
from several 10^ yr to a few Gyr to cover the full ~ 15 Gyr of galaxy ages). The Padua stellar 
evolution library (Bertelli et al. 1994) is used for this purpose. The basic idea is that any SFR can 
be simulated as a linear combination of partial models. 

The practical way to derive the SFR is as follows. A number of boxes are defined in the 
CMDs (observed and partial models) and the number of stars are counted inside them. The 
boxes are defined in such a way they sample stellar evolutionary phases providing information 
about different ages. Let us call N° the number of stars in box j in the observed CMD and NJ^ 
the number of stars in box j of partial model i (the partial model covering the ith age interval). 
Using this, the distribution of stars in boxes of an arbitrary SFR can be obtained from a linear 
combination of the NJl, by 

N^ = kY,a,N^,. (9) 

i 

The corresponding SFR can be written as 

V'(t) =A:^ai^pAi(t), (10) 

i 

where are the linear combination coefficients. A; is a scaling constant, and Aj(t) = 1 if t is inside 
the interval corresponding to partial model i, and Ai(t) = otherwise. Finally, the ip{t) having 
the best compatibility with the data can be obtained by a least-squares fitting of NJ^ to N^, the 
Ui coefficients being the free parameters. 

The finite time interval of partial models introduces a limitation on the time resolution of 
the SFH. But it must be noted that the final limitation is indeed imposed by the accuracy of 
observations and the precision of stellar evolutionary models and also that increasing the time 
resolution of the computed solution is only a matter of computer time. 

Ten partial models have been considered for the analysis of Draco SFH, covering the age 
interval from 15 to 2 Gyr. The latter has been assumed the youngest age present in the CMD 
because it is the age of the isochrone having the turn-off point at "V" ~ 1.5, corresponding to 
the upper point of the Draco blue plume. On the other hand, 16 boxes have been defined in the 
CMD, as shown in Fig. to characterize the distribution of stars. Boxes 1 to 6 sample the 
younger population, living in the upper MS. The consequences of this assumption are discussed 
below. In this analysis we will assume that all these stars are MS stars and not BS. Boxes 7 to 10 
sample the subgiant (SG) region. The information provided by these stars is complementary to 
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the former one. Fainter stars are not considered because they are affected by bigger observational 
effects. In fact, this will limit the time resolution for the oldest stellar ages. Boxes 11 and 12 
sample the HB. These boxes sample the old stellar population. In particular, if age is the second 
parameter, box 11 would contain the oldest pure population. Here, we will assume that this is the 
case, but it must be kept in mind that the nature of the second parameter is not fully confirmed 
and that the blueward extension of the HB depends on stellar evolution parameters that are not 
well controlled. The color uncertainty produced by this is the reason why we define boxes 11 and 
12 to be quite extended in color. 

Box 13 samples intermediate-age to young stars in the core He-burning phase. They should 
have a counterpart in the blue plume and the upper SG region. Region 14 samples the RGB, 
populated by old and intermediate-age stars born in fact in the full, 2-15 Gyr interval considered 
here, as well as some AGBs in the same age interval. Although it provides nothing in terms of 
age resolution when compared with the information contained in other boxes, the number of stars 
in box 14 is a strong tool to normalize the full SFR: independently of the age distribution, the 
integrated SFR for the old and intermediate-age stars must be compatible with the number of 
stars populating box 14. Box 15 samples intermediate-age AGBs. It should have a counterpart in 
the blue plume and the SG region. Finally, box 16 sample young AGBs. The fact that it contain 
no stars indicates that the SFR in the last several million years has been zero or negligible. Linear 
combinations of partial models predicting stars in this region will be automatically rejected. 

We have so far commented qualitatively on the different age intervals sampled by each box 
defined in the CMD. However, the solution for the SFR is found in a global way, considering the 
number of stars in all the boxes. In practice, it is not reasonable to look for the best solution 
(the linear coefficients best reproducing the distribution of stars in the boxes of the observational 
CMD), but to search all the solutions providing a stellar distribution compatible with the observed 
one within some interval. Here, we have adopted as good solutions all those producing a number of 
stars in each box within 2a from the observed value, a for each box being the maximum between 
(iVp^i/^) and 1. The reason for adopting a minimum value of 1 for a is to dump small-number 
effects on the solution. 



The resulting SFR is shown in Figure |lj. Error bars represent the sigma value of all the 
linear combinations producing good solutions for each age interval. The most relevant fact of 
the solution is that the dominant population is old (older than 10 Gyr). It is also remarkable 
that the SFR extends down to a few Gyr-old stars, including a jump between 2 and 3 Gyr. The 
later depends on the assumption we have made that the blue plume is made up of normal MS 
stars. If this is wrong, the SFR of intermediate- age to young stars would be lower than the values 
obtained. However, the intermediate-age star formation activity is also supported by the presence 
of subgiants and core He-burning stars in the CMD. 



This solution cannot be however considered as completely deterministic. In fact, error bars 
are in general comparable to the values of ij^it) for ages between 3 and 10 Gyr. This, however. 
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does not mean that we know very little about the solution, but that we cannot determine with a 
high degree of confidence whether the most conspicuous population is 12-13.5 or 13.5-15 Gyr old, 
for example, or whether the jump at 2-3 Gyr is as strong as it appears in Fig. ^ or is more (or 
less) diluted in time. 

Summarizing, the fact that the dominant population in Draco is old may be considered as a 
well established result of the present analysis. Also, under the assumption that the stars populating 
boxes 1 to 6 (Fig. |l^ are MS stars rather than blue stragglers, a low rate, intermediate-age stellar 
population exist. The jump at 2-3 Gyr is compatible with an overpopulation of stars in box 4 
but mainly comes from the stars populating box 13. These stars can well be ~ 1.2-1.5 Mq in the 
phase of core He burning. 

The SFR obtained here has been used to calculate the total mass in stars and stellar remnants 
(M*) for the inner T.b and 30' (semi-major axis). For this, it has been assumed that a fraction 0.8 
of the integral of the SFR remains locked into stars or stellar remnants. The value obtained for 30' 
can be considered as a good estimate of the total real value, since the fraction of stars at distances 
larger than 30' is negligible (see Fig. ^). Both values are quoted in Table ^ and the second is used 
to calculate the dark matter fraction, k. 

5.2. The subgiant method 

The turn-off is the most sensitive indicator of stellar age for single stellar populations. 
However, in systems such as galaxies, where the stellar population is composite, the turn-offs of 
successive stellar generations mix with each other and with the main sequence of younger stars. 
In such cases, the subgiant (SG) region is preferable for measuring stellar ages. It also has the 
advantage of being, for a given age and metallicity, above the corresponding turn-off point, thereby 
allowing a better signal to noise in the magnitudes. 

Mighell (1997) introduced a simple method based on star counts in the SG region to determine 
the SFR. For this purpose, the SG region is divided with isochrones of a given metallicity and 
increasing ages. Then the number of stars are counted between each two isochrones, which provides 
a present-day age distribution. Normalized by the stellar lifetime in the SG region as a function 
of age, this distribution results in the SFR. In the present case, we perform this normalization 
dividing the observed, present-day age distribution by the distribution corresponding to a synthetic 
CMD computed for a constant SFR and the metallicity adopted for Draco. 

Resolution is worse for older ages, due both to the increasing photometry errors and to the 
decreasing width of the magnitude interval between isochrones. The method provides a good time 
resolution up to ~ 8 Gyr. 

A severe limitation of the method, as it is presently applied, is that it allows only a constant 
metallicity for the stellar population. This, however, is not very important for very low metallicity 
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galaxies such as Draco. Also, as a matter of fact, this method uses only a part of the information 
that is used by the partial model method, for which the solution should be considered less accurate 
and conclusive. However, in the present case, in which metallicity evolution is not tested, the 
subgiant method provides an alternative approach in which the solution rest only on the SG stars. 

In the present analysis, and for the same reasons as given in §4.1, the adopted metallicity 
was Z = 0.0005, with an He abundance of 1" = 0.25. With this input and the Kroupa et al. 
(1993) IMF, a synthetic CMD has been computed for constant SFR in the age interval 1-15 
Gyr. Observational effects have been simulated as explained in Aparicio & Gallart (1995). This 
is important because errors and completeness significantly change along the SG region. The 
simulation has been done six times with different random number seeds. This allows calculating a 
dispersion for the synthetic star counts in each age interval. 

The SG region of Draco's CMD is shown in Figure Isochrones interpolated from the 



Padua library (Bertelli et al. 1994) are shown. The resulting SFR is presented in Fig. IC, which, 
together with a dominant old population, shows an intermediate-age population as far as 2 Gyr 
ago. To estimate the errors two components have been considered. The first is just the square 
root of the star counts for each age interval. It would account for the fact that we work with a 
limited number of stars. The second is the dispersion found for the synthetic star counts and 
would account for the random fluctuations in the observational effect simulation. 

The most interesting thing is that both the partial model and the subgiant methods produce 
qualitatively similar results, including the young bump which, although less clear, is also visible in 



Fig. |16|, slightly shifted to older ages. This stress the stability of the solution and reinforces the 
capability of both methods to obtain the SFH of nearby galaxies. The main difference is that the 
SFR derived by the subgiant method is more softly decreasing and it maintains a higher value 
at intermediate ages. This could be due to observational effects that would spread the stars in 
the subgiant region, the overall result being that fainter stars are moved upwards. In any case, 
the fact that the partial model method uses more information from the CMD eventually makes 
preferable the result obtained with it. 



6. Conclusions 

We have presented a photometric and morphological analysis of the Draco dwarf spheroidal 
galaxy, as well as a study of its stellar content and star formation history based on wide field 
photometry covering about 1° x 1°. 

Two methods [partial model and subgiant) have been used to investigate the SFH of Draco. 
The most important result is that the stellar population is fundamentally old (older than ~ 10 
Gyr). Some 90% or 75% (depending on the method) of the star formation in Draco took place 
before ~ 10 Gyr ago. After that, the star formation continued at a low rate and, apparently, a 
small burst was produced some 2-3 Gyr ago. No star formation activity is detected in the analysis 
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of the CMD for more recent epochs. 

Several morphological and global parameters have been calculated for Draco. Its distance is 
d = 80 lb 7 kpc and its metallicity, [Fe/H] = —1.8 it 0.2. No metallicity gradient is detected along 
the galaxy. 

The star surface density distribution of Draco can be fitted by a single exponential law of 
scale- length a = 5'.0 =b Ofl, for galacto-centric distances larger than 5'. Within errors, the same 
scale is found for the old (older than ^ 9 Gyr) and the young (~ 2-3 Gyr) stars, indicating 
either that both popTilations were formed under the same kinematic conditions or that any initial 
difference was afterwards erased. 

Evidences have been found of the galaxy extending as far as 42' (1 kpc) from its center. 
This is in good agreement with the discovery by Piatek et al. (2001) of stars belonging to Draco 
beyond the previous tidal radius estimate by Irwin k, Hatzidimitriou (1995) (28'). However, the 
absence of plateaus in the density profiles at large galacto-centric distances and the fact that the 
scale-length along the semi-major axis and semi-minor axis (re-scaled to semi-major axis units 
using the cUipticity of the galaxy) are the same within errors, indicates that a tidal-tail is not 
present. More data are probably required to reach stronger conclusions in this sense. 

The central surface magnitude is estimated to be jjniy,, = 24.4 it 0.5 and the core radius 
rc = 7'.5 ± 0'.3 or = 175 ± 7 kpc. 

The possibility suggested by Kroupa (1997) and Kessel &; Kroupa (1998) that the high 
velocity dispersions of dwarf spheroidal galaxies could be accounted for by adequate spatial 
orientations rather that being produced by large amounts of dark matter has been checked. The 
conclusion is that this cannot be the case for Draco. An upper limit to Draco's size along the line 
of sight is ~ 14 kpc; less than 1/3 the value required by Kroupa (1997) models. 

We are grateful to Dr. Kroupa, who made the results of their tridimensional models available 
to us and for several fruitful discussions and to Dr. Gomez-Flechoso, for several useful sugestions 
and discussion on the dynamical part of the paper. 
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on the island of La Palma by the ING in the Spanish Observatorio del Roque de Los Muchachos 
of the Institute dc Astrofi'sica de Canarias. This research has made use of the NASA/IPAC 
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Table 2. Globular clusters used for the metallicity calibration in the 

[Mr, {B - R)] plane 



Cluster 


Metallicity 

ZW 


CG 


Distance 


E{B - V) 


NGC 1904 


-1.7 


-1.4 


15.59 


0.01 


NGC 4590 


-2.1 


-2.0 


15.19 


0.05 


NGC 6723 


-1.1 


-1.0 


14.85 


0.05 



Table 3. Coefficients of the metallicity scales calibration 



Metallicity scale 


a 


b 


c 


ZW 


1.924 


-5.199 




CG 


-3.013 


13.192 


-15.429 
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Tablc 4. Density profile scales 



Stellar population 


«(') 


All stars, full annuli 


4.7 ±0.3 


All stars, semi-major axis 


5.3 ±0.5 


All stars, semi-minor axis^"*"^ 


5.2 ±0.4 


MS, full annuli 


5.0 ±0.1 


MS, semi-major axis 


4.8 ±0.3 


MS, semi-minor axis^^^ 


5.2 ±0.4 


RGB, full annuli 


5.9 ±0.4 


RGB, semi-major axis 


5.1 ±0.6 


RGB, semi-minor axis^^^ 


5.9 ±0.5 



Note. — Parameters given for semi-minor axis have been 
scaled onto semi-major axis units; i.e., they have been divided 
by h/a. 



Table 5. Morphological and integrated parameters of Draco 



Parameter 


Value 


"00 


1 Thonm o 
11 2Xj .6 


t>00 


57 55 


H"V",o 


0/1 A \ n ^ 
Z4:A dz (J.i3 




nine: 




26.9 ± U.O 


a 


o'U ± U:l 




7.5 ± U.O 




(175 ± 7 kpc) 


d 


8U ± 7 kpc 


[Fe/H] 


-1.8 ±0.2 




10.1 ±0.5 




10.6 ± 0.5 


Rt 


9.6 ±0.5 


Miiyif 


-9.4 ±0.5 


Mb 


-8.9 ±0.5 


Mr 


-9.9 ±0.5 


L"v" (Lq) 


4.7 X 10^ 


MvT {Mq) 


22 X 10^ 


Mvt/L^^V" {Mq/Lq) 


47 


M,ia < 7'.5) (Mq) 


2.1 X 10^ 


M^a < 30') (Mq) 


4.3 X 10^ 


M^{a < 42') (Mq) 


5.3 X 10^ 


K=l- {M^{a < A2')/Mvt) 


0.98 
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Table 6. Comparison of morphological parameters 



Reference 


rcC) 


nC) 


Hodge (1964) 


6.5 


26 ±2 


Irwin & Hatzidimitriou (1990) 


9.0 ±0.7 


28.3 ± 2.4 


Lake (1990) 


7.0 ±2.5 


55 


Pryor & Kormendy (1990) 


7.1 


36 


This work 


7.5 ±0.3 


42 
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Fig. 1. — Digitized Sky Survey image of Draco. The observed fields are shown. North is up, East 
is left. 
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Fig. 2. — Isopleth map of Draco taken from Irwin k, Hatzidimitriou (1995). Coordinates refer to 
the scales given by these authors. The observed fields are shown. North is up, East is left. 



-28- 



1 - 



0.8 - 



0.6 - 



0.4 - 



0.2 - 







— I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — r 

■ i i i 



J I L 



J I I I I I I L 



20 21 



22 23 
"V" 



J I I L 



24 25 



Fig. 3. — The completeness factor as a function of "V" for field 1. 
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Fig. 4. — The CMD of Draco. Stars within the inner 30' (semi-major axis) are plotted. Beside a 
well populated MS and a steep RGB, indicative of relatively low metallicity, the most significant 
features are the extended blue plume and the mainly red HB. Some blue HB stars and an AGB 
sequence are also visible. 
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Fig. 5. — CMDs of four elliptical annuli centered on the center of Draco and of semi-major axis as 
indicated. 
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Fig. 6. — The upper part of the Draco RGB. Only stars in the inner 10' have been plotted. The 
thick track shows the line of constant absolute / magnitude for M/ = —3. The curve has been 
calculated using a synthetic CMD. The scale plotted on the line shows the values of {V — I)/=_3 
along it. The fiducial RGB of Draco is shown by the thin line. The point where this line crosses 
the curve provides {V — I)i=-3 for Draco, from which the metallicity can be estimated. 
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Fig. 7. — The fiducial RGBs of the three globular clusters used to set the metallicity calibration 
in the [Mr, {B — R)] plane discussed in the text are shown superimposed on the CMD of Draco 
(inner 10'). This shows that the metallicity of Draco is probably between those of NGC 4590 and 
NGC 1904. 
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Fig. 8. — (B — R)r=-2A5 (the color of the RGB measured at = —2.45) of Draco stars for 
elhptical annuh of increasing size. If a metalhcity gradient should exist in Draco it would be put 
into evidence by a gradient of {B — R)r=-2A5- 
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Fig. 9. — Surface star density and surface magnitude distributions in Draco. All the resolved stars 
have been considered to plot the upper panel figure, while only MS and RGB stars have been 
respectively used for the medium and lower panel figures. Stars have been counted in elliptical 
annuli of eccentricity 1 — a/b = 0.29, position angle 82° and increasing semi-major axis. Horizontal 
axis gives the semimajor-axis in arcmin. Left-hand vertical axis gives the logarithm of the star 
density per arcsec^. The right-hand vertical axis gives an estimate of the surface magnitude, 
calibrated as described in text. Open circles represent the raw density distribution. Filled circles 
show the density distribution after subtraction of the sky foreground contamination, estimated as 
described in text. 
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Fig. 10. — Surface star density of stars along the semi-major and semi-minor axis of Draco, after 
correction of foreground contamination. MS (upper panel) and RGB (lower panel) stars have been 

considered. The density distribution obtained for star counts in the full elliptical annuli are also 
plotted for comparison. Straight lines show the fits done using full elliptical annuli distributions. 
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Fig. 11. — The spatial distribution of stars according to RSl-4 model of Kroupa (1997), which 
reproduces the velocity dispersion of Draco. The Z axis corresponds to the line of sight. Coordinates 
are given in kpc and are centered on the center of Draco. 
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Fig. 12. — The left-hand panel shows a synthetic CMD for an old stellar population, similar to 
the main population found in Draco (see §4). The right-hand panel shows the same population 
including the simulation of the distance dispersion corresponding the RSl-4 model of Kroupa (1997), 
as shown in Figure 
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Fig. 13. — The boxes defined to sample the stellar population of Draco, overplotted on the CMD. 
These boxes are used to derive Draco's SFR with the partial model method (see text). 
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Fig. 14. — The solution for the SFR of Draco using the partial m,odel method. Two sohitions are 
given. The first for the inner 7'.b and the second for the region between 7'.b and 30'. Error bars 
account for the dispersion of possible solutions for each age range. 
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Fig. 15. — Subgiant region of Draco's CMD. Isochrones of several ages are overplotted, 
corresponding to the ages used to derive the SFR with the subgiant method (see text). From 
up to down, ages are 1.25, 2, 3, 4, 5, 6, 8, 10, 12 and 15 Gyr. Diagonal straight hues show the color 
interval used to limit the SG region. 
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Fig. 16. — The solution for Draco's SFR using the subgiant method. The inner 10' (semi-major 
axis) have been considered. 



